Transport coefficients of a unitarized pion gas 



Juan M. Torres-Rincon 

Dept. Ffsica Teorica I, Universidad Complutense de Madrid, Spain 

November 17, 2011 

(N ' 

> 
O 

The latest experimental results in relativistic heavy-ion collisions show that the matter there 
produced requires transport coefficients because of the important collective properties found. We 
review the theoretical calculation of these transport coefficients in the hadron side at low temper- 
atures by computing them in a gas composed of low energy pions. The interaction of these pions 
is taken from an effective chiral theory and further requiring scattering unitarity. The propagation 
of D and D* mesons in the thermalized pion gas is also studied in order to extract the heavy quark 
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Abstract 



r~| | diffusion coefficients in the system. 
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1 Introduction: Transport Coefficients 

t^. . Transport coefficients are quantities that describe how the system tends to equilibrium when an external 
^} ', perturbation has been applied to it. This perturbation produces a gradient of some hydrodynamical 
field (pressure, temperature, velocity...) followed by a response flow in the system (Le Chatelier princi- 
ple). This response flow pushes the system towards the equilibrium state. At first order in gradients, 
the proportionality coefficient between the hydrodynamical gradient and the response flow is called a 
transport coefficient (TC): 

>< 

5_i ■ TC x V (Hydrodynamics Field) = —Response Flow . (1) 

d | 

The response flow corresponds to the motion of a conserved quantity density so there exists one 
transport coefficient related to each conserved quantity of the system. Some classical examples of 
transport coefficients are the shear and bulk viscosities rj, ( both related to the momentum conservation, 
the heat or thermal conductivity k related to the energy and particle conservation and the electrical 
conductivity a if an electric charge is present in the system. 

In a relativistic heavy-ion collision, some of them can be extracted from those observables encoding 
the collective behaviour of the plasma. One of the most remarkable transport coefficient extracted in 
this way is the shear viscosity, usually normalized to the entropy density r]/s. 

2 Flow Harmonics and Viscosities 

The way to estimate the value of the viscosity over the entropy density is to compare the experimental 
measurement of a collective observable to the result given by hydrodynamics simulations (for which 
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Figure V. Left panel: Differential elliptic flow as a function of the transverse momentum as measured 
at RHIC (dots) and extracted from the hydrodynamic calculation in [2] (lines) for several values of r]/s 
and Glauber initial conditions. Right panel: pr-integrated elliptic flow as a function of the centrality 
as measured at RHIC and LHC compared to the extracted v 2 from the hydrodynamic code of [3] for 
the optimal values of rj/s and (/s and Glauber initial conditions. Figures courtesy of M. Luzum and 
P. Romatschke (Copyright 2008 by The American Physical Society) and P. Bozek, respectively. 



r]/s is an input parameter). The most impressive collective phenomenon is the presence of momentum 
anisotropies in the expanding fireball due to initial spatial anisotropies in noncentral collisions. The 
quantitative effect of these anisotropies is reflected in the flow harmonics v n , defined as the averages 
over the measured particles in some event: 

v n = (cosn(0i - *hp)) , (2) 

where <pi is the azimuthal angle of the outgoing particles and ^rp is the reaction plane, defined by the 
beam axis and the impact parameter vector between the two incoming nuclei. For a recent estimation 
of the flow harmonics made by the ALICE collaboration see pp. 

For noncentral collisions, the most prominent flow harmonic is the elliptic flow V2 and it is quite 
dependent on the value of r]/s. By matching the experimental behaviour of t>2 and the hydrodynamic 
calculations one can estimate the approximate value of rj/s. This value turns out to be very close to 
77/s ~ l/(47r) (see left panel of FigJT]). The value of the bulk viscosity over entropy density {(/s) can 
also be estimated following the same method (see right panel of Fig. [p. 

The newest investigations along these lines are obtaining more precise values for the two viscosities 
in order to restrict other hydrodynamic parameters and discriminate among different initial condition 
models. As an example, the more detailed calculations have began to include the entire temperature 
dependence of the shear viscosity (see for instance [I]) and of the bulk viscosity (see for example [5]). 
We are going to describe the theoretical calculation of these transport coefficients and its temperature 
dependence in the low temperature phase, i.e. the hadronic medium. 

3 Hadronized Medium: The Pion Gas 

We will consider that the matter produced after the hadronization time is mainly composed by pions. 
Experimentally, it has been found that the pion multiplicity is of the order of ten times larger than 
the total multiplicity of the next abundant degree of freedom, the kaon. At moderate temperatures 
T ~ 140 MeV the next mesonic states (kaons and 7] mesons) should be included in the calculation. 

The properties of a pion gas at low temperatures are not calculable from perturbative QCD. We 
will use effective field theory (ChPT) to describe the interaction between pions. The SU(2) chiral 
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perturbation theory (ChPT) provides a good description of the scattering amplitudes in the energy 
range we are interested in. 

We will assume that at low temperature the pion gas is dilute enough to only consider binary elastic 
collisions. The suppression of the number changing processes (2 — > 4 or 4 — > 2) comes from two sources. 
First, the Lagrangian terms describing these processes are of higher order in momentum and therefore 
are polynomially suppresed at low temperature (as long as p ~ \frnT, where m is the physical pion mass 
m = 138 MeV). Second, in the 2 — > 4 process the energy needed to create two more pions makes this 
scattering suppressed by a Boltzmann factor of e~ 2m . Finally, the process 4 — > 2 is evidently opposed 
to the spirit of Boltzmann assumptions for molecular chaos and does not occur at low density. The 
encounter of four pions is manifiestly unlike to occurs in a dilute gas. Including only elastic scattering, 
the pion number is an effectively conserved quantity and a (pseudo-) chemical potential for the pion 
must be introduced. 

The Lagrangian of ChPT with physical pion masses follows a systematic counting based on terms 
with even powers of pion momentum and pion mass over the scale A = AnF ~ 1 GeV: 

C = £ 2 + U + - • • , (3) 

where £2™ is 0[(p/A) 2n ] and 0[(m/A) 2n }. The G— parity conservation guarantees that 2n is an even 
number and it amounts to have interactions between an even number of pions. 

The lowest order contribution to the scattering amplitude is the elastic scattering of two pions. We 
want to stress again that inelastic channels are suppresed at Lagrangian level (6-particle interaction is 
described by C§) and because of Boltzmann suppression in the final phase-space. 

From the power counting scheme in the ChPT Lagrangian the partial scattering amplitudes tjj 
(labelled by the isospin and angular momentum channels) are essentially polymonials in momentum: 

tff - o( P 2n ) , (4) 

where this term comes from the part of the Lagrangian of order £2™- Truncating the expansion at some 
finite term gives a scattering amplitude that increases with energy as tfj^ ~ s n . This effect causes an 
unphysical increase of the pion-pion cross section at moderate energies and leads to a breaking of the 
exact unitarity condition for the scattering amplitude. This conditions reads: 

Im tjj(s) = sfl - Am 2 /s \tjj(s)\ 2 . (5) 

The partial amplitudes in (T4]) only satisfy this equation perturbatively, e.g. Im tfj(s) = 0, Im £/j(s) = 
1 -Am 2 js \tf)(s)\ 2 ... 



We use the Inverse Amplitude Method [6],j7], a dispersive method that constructs a new partial 
amplitude tu(s) from the pertubative ones that exactly satisfies the unitarity condition 03). At lowest 
order, the unitarized amplitude turns out to be: 

This new amplitude, as rational combination of perturbative amplitudes, makes the pion-pion cross 
section saturate at moderate energies and allows to dynamically generate the p and a resonances in the 
IJ = 11 and IJ = 00 channels, respectively. 

For calculating the transport coefficients, one needs to solve a transport equation for the one-particle 
distribution function f w . This transport equation is of the well-known form: 

-^ = c[r,r\. (7) 
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For a gas of pions the explicit form of the right-hand side is that of the Boltzmann-Uehling-Uhlenbeck 
(BUU) equation: 



dJl = g JL 
dt 2 



/ dr 12t3p [fif-{\ + ft){i + /;) - / 3 v;(i + /na + m 

where g n = 3 and the measure c£r 12i 3 P constains the details of the pion interaction (1, 2 — > 3,p) 
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dTl2 ^ =2E p lTl2 P (^j^ 27 ^^ 1 + k 2-h- P ). (9) 

The integral equation ([8]) is linearized by using the Chapman-Enskog procedure, consisting in an expan- 
sion of the distribution function in powers of the Knudsen number (e = mean free path / characteristic 
size of the inhomogeneities of the system) 

f; = fBose + ef^ + --- (10) 

In this sense, the perturbative solution is only valid if hydrodynamics can be applied to our system. 
Once the transport equation has been linearized an approximate solution for is obtained expanding 
this function in powers of the variable E p /m, 

N / TP \ n 

/ (1) = E^p) , (ii) 

up to some finite N. Finally, the linearized BUU equation must be numerically inverted to extract the 
coefficients A n . 

Once is built from the A n coefficients, the transport coefficient is expressed as a simple integral 
(or a scalar product): 

TC = ([ V(Hydrodynamic field) • ( 12 ) 

The shape of the hydrodynamic fields and the parametrization of are dependent on the transport 
coefficient we are calculating. The details of the calculation for the shear viscosity and bulk viscosity 
can be found in |8]-[9] and [TO], respectively. 

In Fig. [2] we show the results for the transport coefficients of a pion gas as a function of temperature 
and for several values of the pion chemical potential (/i < m). In the top panel we show the shear 
viscosity, the shear viscosity over entropy density and the bulk viscosity. In the bottom panel we show 
the bulk viscosity over entropy density, the thermal conductivity and the electrical conductivity. 



4 Adding Heavy Quarks to the Medium 

We now include a heavy degree of freedom (c quarks for concreteness) inside the thermalized pion gas. 
The charm quark -hadronized into a D or a D* meson- propagates in the hot medium and interacts with 
it transferring energy to the pions due to these collisions. The charm diffusion coefficients describes the 
energy loss and momentum spreading of these heavy mesons. 

Experimentally, the diffusion coefficients are important in order to describe the nuclear modification 
factor, Raa = ^^ydN^/dpT anc ^ ^ e e ^ n P^ c ^ ow °f ^ ne electrons coming from the semileptonic decays 
of these heavy mesons [UJ. For more details we refer the reader to [13J. 

The interaction of the D and D* mesons with the pions is described by an effective Lagrangian. This 
Lagrangian encodes the chiral and the heavy quark symmetries because of the presence of the pions 
and the charmed mesons, respectively. These two symmetries are the guiding principles to construct 
the sucessive orders of the effective Lagrangian [T2].[T3]. 
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Figure 2: Transport coefficients of a unitarized pion gas. Top panel: Shear viscosity, shear viscosity over 
entropy density and scaled bulk viscosity. Bottom panel: Bulk viscosity over entropy density, thermal 
conductivity and electrical conductivity. 



Most of the terms in the Lagrangian are preceded by low energy constants that are not known. 
Some of these constants can be fixed by symmetry arguments and others should be adjusted by asking 
them to correctly reproduce the experimental decay widths or cross sections of charmed resonances. 
We work at NLO in the chiral expansion and LO in the heavy quark expansion. After constructing the 
perturbative scattering amplitude for the D — n and D* — n scatterings, we unitarize the amplitudes in 
the same spirit as for the pion gas. 

The unitarization is done by using the "on-shell' unitarization technique [H],[IH] that consists in 
solving the Bethe-Salpeter equation for the D — n rescattering. The final expression of the unitarized 
amplitude resembles very much the rational structure of (j6]). Like the Inverse Amplitude Method, the 
"on-shell" unitarization provides a well-behaved cross section. Moreover, we can dynamically generate 
the .Do (2400) and _Di(2430) resonances, that are automatically incorpored in the calculation. 

The transport equation for the distribution function of the c quark reads 

df c 

^=C[/ C ,r], (13) 

in complete analogy with the BUU equation for a gas of pions ((7|). However the mass of the charmed 
mesons (Mo, Mp*) -much larger than the pion mass (m) or the temperature (T)- provides the scale 
hierarchy 

Mo, Mr>* Cm~T~ transfered momentum , (14) 

that simplifies the transport equation for solving the heavy quark distribution function. The BUU 
equation tranforms into a Fokker-Planck equation under Eq. (TTJ 



df c (t,p) _ d 
dt dpi 



^(p)/ c (t,p) + T^(r, J (p)r(t, P )) 



(15) 
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Under the assumptions of homogeneity and isotropy one finds three different coefficients that depends 
on the heavy meson momentum, one drag force F(p) and two diffusion coefficients T Q (p) and Ti(p): 



Fi(p) = F(p) Pi , (16) 

r.(p) = r (p)^-^)+r l( p)^. (17) 

Because of the fluctuation-dissipation theorem there are only two independent coefficients. The 
explicit expression for these coefficients in terms of the scattering amplitude can be found in [13] . 

In the static limit p — > only one of them is independent because the two diffusion coefficients 
become degenerate and also because the Einstein formula relates this one to the drag force T(p — > 0) = 
F(p -> 0)TM D . 

There are other related quantities, with direct physical interpretation. On the one hand, the spatial 
diffusion coefficient D x -that appears in the description of the Brownian motion- is related to the static 
diffusion coefficient. On the other hand, momentum and energy losses are directly related to the drag 
force. The explicit equations for them are the following: 

D X = T 2 /T; -^ = ExF(p); — = pxF(p). (18) 

We summarize the results in Fig. |3j We show the two diffusion coefficients T (p) and Ti(p) as a 
function of the temperature for selected heavy meson momenta. These diffusion coefficients enter in 
the description of the momentum anisotropies, i.e. the elliptic flow of the heavy mesons. We also show 
the same dependence of the heavy meson drag force, that tells us the momentum losses of the heavy 
meson affecting to the nuclear modification factor of D and D* mesons. In the same figure, we show 
the energy and momentum losses of the heavy meson that for instance, at p = 1 GeV leaves about 50 
MeV of momentum per Fermi travelled. Finally, we show the spatial diffusion coefficient D x together 
with some other results along the same lines and the calculation for a perturbative heavy quark inside 
the quark-gluon plasma (see [13] and references therein). 



5 Conclusions 

We have described the theoretical calculation of the transport coefficient of a meson gas at low temper- 
atures. We have studied the pure pion gas and the case where charm degrees of freedom are added to 
this gas. For both cases, we have used an effective Lagrangian to describe the interactions between the 
mesons and we have provided unitarization schemes to ensure that the scattering amplitudes are prop- 
erly unitarized. All these coefficients provide valuable information about the experimental observables 
that reflects collective phenomena of the medium, such as flow harmonics and nuclear modification 
factors. We believe that a better understanding of these non-equilibrium effects will provide a better 
estimation of these observables and an accurate description of the properties of the matter produced at 
relativistic heavy ion collisions. 

Work supported by grants FPA 2008-00592 and AIC10-D-000582 (Spain). The author is recipient of an 
FPU scholarship from the Spanish Ministry of Education (AP2008-0083). Work done in collaboration 
with A. Dobado, F.J. Llanes- Estrada, L.M. Abreu and D. Cabrera. 
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Figure 3: Top panels: Values of the two diffusion coefficients and the drag force as a function of 
temperature and momentum. The three coefficients are related through the fluctuation-dissipation 
theorem. Bottom panels: Results for the energy loss, momentum loss and the spatial diffusion coefficient 
in comparison with other approaches. All the plots has been extracted from [13] . 
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